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Abstract 



We show how to divide a coupled multi-spin system into an small subset of "active" 
spins that evolve under chemical shift or scalar coupling operators, and a larger subset of 
"spectator" spins which are returned to their initial states, as if their motion had been 
temporarily frozen. This allows us to implement basic one-qubit and two-qubit operations 
from which general operations on iV-qubits can be constructed, suitable for quantum 
computation. The principles are illustrated by experiments on the three coupled protons 
of 2,3-dibromopropanoic acid, but the method is applicable to any spin-1/2 nuclei and to 
systems containing arbitrary numbers of coupled spins. 
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1 Introduction 



In this letter we address the problem of how to use NMR techniques to manipulate subsets 
of spins in a multi-spin system while leaving the states of all the other coupled spins 
unchanged. We have in mind the use of these manipulations to generate pulse sequence 
modules which will be useful for NMR quantum computation Jj], |2[ ||L [5|, ||, [?], @, M. 
This leads to an important difference of emphasis from conventional NMR spectroscopy 
where one can assume that the states of the spins before the manipulation are known, 
and indeed one often arranges these states (for example as z-magnetization) so that the 
manipulation has a simple effect on the spins. One is certainly not concerned about the 
effect that the pulse sequence might have on spins which were prepared in different states. 
Furthermore, one is typically only interested in the effect of the manipulation on a subset 
of the spins and it is irrelevant how the other spins evolve. In contrast, in applications 
to NMR quantum computing, the state of the system before the given pulse sequence 
module is assumed unknown. The module is a logic gate which is designed to carry out a 
specified unitary transformation which can be applied to any input, the remaining qubits 
must remain passive spectators. 

A system of N spins evolves under a Hamiltonian which includes chemical shifts for each 
spin and scalar coupling terms for every pair of spins. Thus if one wishes to manipulate a 
particular spin, and the manipulation is not effectively instantaneous, all the other spins 
will evolve, in a complicated, coupled, manner. We will show here how to apply patterns 
of pulses so that only the desired chemical shifts and scalar couplings affect the evolution 
of the system as if we had "switched off" all the other terms in the Hamiltonian. For 
example, we will show how to produce pulses whose net effect is to allow only one spin to 
precess under its chemical shift: no matter what the input states of the other spins are, 
at the end of the pulse sequence they will be returned to their initial states. 

For concreteness we will use a weakly-coupled three spin (ISR) system to demonstrate 
the theoretical and experimental results. This system has the "free" Hamiltonian (i.e. the 
Hamiltonian in the absence of pulses) of the form 

Ho = ujiI z + ljsSz + u R R z + 2irJ IS I z S z + 2irJ IR I z R z + 2-kJ sr S z R z . (I) 

Using patterns of refocusing pulses we will show how to create unitary operations on the 
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spin system as if many of the terms in the Hamiltonian were zero. We will demonstrate 
theoretically and experimentally how to generate three basic operations: 



• Ol Chemical shift evolution of a single spin 

• 02 Rotation of a single spin by an arbitrary angle about a general axis (this gener- 
alises Ol) 

• 03 Evolution under scalar coupling leading to the two-qubit operation whose prod- 
uct operator representation is exp(i(p2I z S z ) 

Straightforward generalisations of the sequences will produce analogous modules to Ol, 
02 and 03 for systems of any number of spins and these modules are sufficient to generate 



any unitary operation on N qubits. This follows from the important result in [10] that 
any arbitrary unitary transformation of N spins can be built up from two simple building 
blocks: these may be taken to be, for example, an arbitrary rotation of a single spin, and 
the two-qubit operation exp(— i^2I z S z ). The reason for this particular choice of basic 
operations is that they may easily be generated using the natural evolution of the system. 
We note that the two-qubit operation exp(— i^2I z S z ), which has matrix representation 
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is closely related to the controlled-not (CNOT) operation on two qubits. 

The first operation Ol is implemented by the pulse sequence given in Fig 1. The sequence 
of refocussing n pulses on spins R and S creates a unitary evolution of the system as if 
the only non-zero term in the Hamiltonian were ujl z . Thus the effect of the sequence is 
to act on the complete system with the unitary operator e l< ^ Iz . As an example of how the 
generalisation to iV-spins works, consider the case of four spins labelled ISRQ. Evolution 
of the I spin alone can be arranged during a period of 16r by having n pulses on S at 
(8m + 4)T, m = ... 1, on R at (4m + 2)r, m = . . . 3, and on Q at (2m+l)r, m = . . .7. 

Returning to the ISR system, we see that we can convert the phase evolution sequence (a 
rotation about the z axis) into a general excitation pulse 02. The /-spin magnetization 
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may be rotated through an angle <fi about a tilted axis in the xz plane by enclosing the 
unitary operator e lr ^ Iz between two hard pulses applied along the — y and +y axes: 

e -i(B(Iy+Sy+Ry) e i4>h e + il3(Iy+Sy+Ry). ^ 

the hard pulses are effectively instantaneous so that no evolution under TLq takes place. 
This tilts the rotation axis away from +z through an angle (5 in the xz plane. The rotation 
can be made completely general by enclosing sequence (||) between two z-rotation Ol 
modules: 

Rotation now takes place about a tilted radiofrequency field in a vertical plane that sub- 
tends an azimuthal angle 7 with respect to the +x axis. As a specific example we will 
demonstrate experimentally the case <j> = it, (3 = it/A, 7 = 0. 

The final basic operation is the two-qubit operation 03. The sequence of refocussing 
pulses is given in Fig 2. We note that this sequence, together with the previous ones, can 
be used to create the unitary operation whose total effect is to perform a CNOT operation 
between I and S (with the R spin unaffected). Recall that the effect of 03 is, in terms of 
product operators, exp(i(j)2I z S z ). The CNOT operation on two spins may be represented 
by the unitary matrix 



t^CNOT 
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(5) 



The corresponding product operator representation is 

It can be seen that each term in the product (except for the irrelevant overall phase) can 
be created using Ol, 02 or 03. 

Thus the pulse sequences we describe in this letter give us a set of basic modules which 
are sufficient to produce an arbitrary unitary operation on a system of 3 qubits. In 
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practice there may be more straightforward pulses which suggest themselves as NMR 
implementations of the essential parts of given gates. However there are often "corrections" 
which need to be made to implement the exact unitary operator required |2], |9| . Consider 
for example the controlled-controlled-not (CCNOT) gate on three qubits: 
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(7) 



This may be implemented with a sequence of five controlled rotations [1C]. However using 
a line-selective pulse it is straightforward to generate the key part of this gate || , giving 
the unitary operator 
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(8) 



It is not difficult to show that 



^CCNOT = e*4 J * 



i * 2IxS * L^CCNOT, 



(9) 



up to an irrelevant overall phase. Thus the exact unitary operation required, i/cCNOTj 
can be implemented by performing t/ccNOT followed by a sequence of corrections of the 
type we describe in this letter. 



2 Experimental realization 

We take as representative practical example a homonuclear weakly-coupled three-spin 
(ISR) system where all three coupling constants Jis, Jib, and Jrs are non-zero and 
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resolvable. Relaxation effects (spin-spin and spin-lattice relaxation and the nuclear Over- 
hauser effect) are specifically excluded from the discussion since they change intensities 
in an irreversible manner; of course if the pulse sequences have appreciable duration, 
relaxation effects will be observed. Multiplet-selective radiofrequency refocusing pulses 
are employed in a symmetrical arrangement to refocus the appropriate chemical shifts and 
spin-spin splittings. An essential property of these soft pulses is that they should act as 
"universal rotors" , that is to say, within their effective bandwidth, they induce the same 
rotation whatever the initial state of the nuclear magnetization. Gaussian shaped pulses do 
not satisfy this requirement; instead time-symmetric, shaped radiofrequency "RE-BURP" 
pulses [11 1 are used. These refocusing "7r" pulses are never applied to two different nuclear 
spins simultaneously, since that would introduce spin-echo modulation through the scalar 
coupling and create multiple-quantum coherence through the double-resonance two-spin 
effect (TSETSE) Jl2|| . The refocusing pulses are employed in a sequence designed to 
take advantage of the fact that radiofrequency pulse imperfections are corrected on even- 



numbered spin echoes [13|. 

Experiments were carried out at 22C on a Varian 400 MHz high-resolution spectrometer 
on the three coupled protons (R, I and S) in 2,3-dibromopropanoic acid dissolved in 
benzene-dg, where Jis = —10.1 Hz, Jjr = 11.3 Hz, Jrs = 4.3 Hz, 5jr = 188.5 Hz 
and Sis = 219.5 Hz. After removal of dissolved oxygen by bubbling argon through the 
solution, the measured spin-lattice relaxation times were 9.5 s (R) and 1.8 s (J and S). 
Although this particular homonuclear system provides a simple example to demonstrate 
the principles involved, it is not ideal for the purpose. The chemical shift differences are 
rather small, placing heavy demands on the selectivity of the soft ir pulses, the uniformity 
of excitation across a chosen spin multiplet, and the effectiveness of the suppression off- 
resonance. Some strong coupling effects are also evident, particularly between spins R 
and /. Analogous experiments on heteronuclear systems (or on homonuclear carbon-13 
systems) would be expected to pose fewer practical problems. Selectivity would no longer 
be crucial, and because the pulses would be of much shorter duration, relaxation and 
coupling effects during the pulses could be safely neglected. 
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2.1 Module Ol: Selective phase evolution 

The basic module Ol is designed to refocus the spin-spin splittings of the active (I) spin 
but allow chemical shift evolution to accumulate a phase shift (ft, leaving the coupled 
spectators R and S unaffected (Fig. 1). After excitation of magnetization by a hard 
pulse applied about the x-axis, the divergence of /-spin magnetization trajectories due to 
Jir and Jis is refocused by soft ir pulses applied to R and S, so that only the chemical 
shift effect persists. It was found that the effects of pulse imperfections were reduced by 
applying the two S pulses with the same phase (+y) and by setting the phases of the four 



R pulses to (+y + y — y — y) according to the MLEV-4 prescription [14|. The phase 
evolution diagram (Fig. 1) illustrates the refocusing process. The S and the R spins remain 
passive spectators because their chemical shifts and spin-spin splittings are refocused by 
the 7r pulses, and certain pulse imperfections are compensated by the operation of the tt 



pulses acting in pairs [13|. At the end of the sequence the I spins accumulate a phase 
shift (ft = 2tt5i8t radians which can be controlled by adjusting either the duration of the 
sequence (8r) or the offset 5i Hz between the /-spin chemical shift and the transmitter 
frequency. 

There is a minor instrumental complication when soft radiofrequency pulses are used in 



this manner, since each radiofrequency field induces a small Bloch-Siegert shift [15] on 
neighbouring resonances, always in the direction to move the neighbour resonance away 
from the irradiation field. This translates into a small additional increment in the phase 
of the nuclear precession. For the I spins these small phase errors merely accumulate and 
can be corrected by a slight readjustment of the overall duration of the sequence. For 
the S spins the symmetry of the sequence is such that these phase errors are effectively 
refocused by the soft n pulses. Unfortunately the perturbation of the motion of the R spins 
by the pulses applied to S is not refocused. The Bloch-Siegert effect on R was therefore 
compensated by applying "ghost" ir pulses at the same time as the S pulses and with 
the same amplitude and duration, but with the offset from the R resonance reversed in 
sign. These ghost pulses fall in empty regions of the spectrum, exciting no NMR response 
directly. With this precaution an experimental test of the e l?Uz module gave the spectra 
illustrated in Fig. 3. The duration of the sequence (8r) was incremented so that the 
accumulated phase (ft of the /-spin response varied over a range of 2ir radians. Note the 



6 



attenuation of the I-spin signal due to T 2 * losses over the relatively long durations (8r 
varied between 1.15 and 1.23 seconds), whereas this effect is refocused for the R and S 
responses. 

A "do-nothing" module that refocuses all chemical shift evolutions and spin-spin splittings 
can be fashioned from the module which creates e l ^ Iz by choosing 8tSi to be an integral 
number of complete rotations, (as is essentially the case in Fig 3(a) and 3(e)). A general 
"do nothing" sequence, valid for any duration 8r, can be constructed by introducing two 
hard (non-selective) 7r pulses at times 4r and 8r, where the trajectories of /, R and S are 
at a focus (Fig. 1). 

2.2 Module 02: Selective rotation 

The Ol phase evolution sequence (a rotation about the z axis) can be converted into a 
module which produces rotation of a given spin by an arbitrary angle about an arbitrary 
axis (with no effect on the other spins). The principle of "decoupling" the rotation of the 
active (I) spins from the motion of the spectator spins (R and S) was tested on the protons 
of 2,3-dibromopropanoic acid by imposing an /-spin rotation through <f> = it radians about 
an axis in the xz plane tilted at /3 = 7r/4 radians. This takes +z magnetization to the +x 
axis, for example. These experiments exploited an improved RE-BURP pulse (effective 
bandwidth 25 Hz), re-optimized to give more uniform re-focusing and better out-of-band 
suppression. Precautions were taken to compensate the effects of spin-lattice relaxation 
(which acts unequally during the different r intervals) by increasing the tilt angle (5 be- 
yond 7r/4 radians. The experimental results (Fig. 4a) confirm the principle of decoupled 
rotation, showing a strong excitation of the /-spin response and residual responses from 
the spectators R and S reduced by more than an order of magnitude. Simulations were 
performed for this sequence based on numerical integration of the Liouville-von Neumann 
equation, neglecting relaxation but allowing for strong coupling. They indicate similar 
levels of residual responses for R and S (Fig. 4b). 

2.3 Module 03: Selective evolution under the scalar coupling 

The sequence shown in Fig. 2 implements the principal part of a CNOT gate, namely the 
unitary operator e l< ^ 2IzSz . Hard 7r pulses are applied at times 2r and 6r, where the Jjr 
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splittings of the /-spin response and the Jsr splitting of the /S-spin response are at a focus. 
In this experiment the preparation stage involved the excitation of y magnetization. The 
/-spin and 5-spin chemical shifts are refocused but the splitting Jis continues to generate 
diverging trajectories throughout the sequence, giving rise to a phase evolution of ±8tttJis 
radians on both the /-spin and 5-spin responses. This is illustrated in Fig. 5 for phase 
increments of ±7r/2 radians. On traces (b) and (d) the (antiphase) I and S responses 
would have appeared in dispersion; for the sake of clarity the receiver phase was changed 
by 7r/2 to reset them to the absorption mode. The R spin response remains essentially 
unchanged throughout. 



3 "Much Ado About Nothing?" 

Whereas classic NMR experiments commonly ignore collateral perturbation of coupled 
neighbour spins, a quantum computer requires that all these "spectator" spins return 
to their initial states at the end of the manipulation. For pulse sequences that are not 
instantaneous, this "do nothing" operation is by no means trivial; it necessitates exact 
refocusing of chemical shifts and spin-spin splittings. 

We have demonstrated how selectivity can be achieved in a homonuclear coupled three- 
spin system for three important basic operations - phase evolution, an arbitrary rotation, 
and evolution under the scalar coupling. Note that these operations do not assume that 
the spin system is initially at Boltzmann equilibrium; the new "modules" are unitary 
operators that act on an arbitrary initial state of the spin system. 

For NMR to be a useful technology for quantum computation, it will be necessary to ex- 
tend current ideas to systems of many more spins. Indeed recent work suggests that only 



for such extended systems will the computations be fully quantum mechanical [16|. The 
extension of these soft pulse sequences to accommodate more spectator spins is straightfor- 
ward, although in practice this lengthens the duration of the sequence. We must remember, 
however, that future implementations of NMR quantum computers need not be restricted 
to groups of coupled protons, so the selectivity requirements could be far less severe than 
in the present proton experiments, implying much shorter radiofrequency pulses. 
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Figure Captions 



Fig. 1. (Top) The Ol module. This sequence of soft it pulses (ellipses) is designed to 
allow phase evolution of the four /-spin magnetization components due to their chemical 
shift (dotted line) without spin-spin splittings, returning the coupled R and S spins to 
their initial states. (Bottom) Phase evolution diagram for the four / lines. 

Fig. 2. (Top) The 03 module designed to focus the /, S and R chemical shifts but 
allow divergence of / and S magnetization components due to Jjs- Ellipses represent soft 
pulses. (Bottom) The phase evolution of the four / lines. 

Fig. 3. Experimental test of the Ol module which allows phase evolution of transverse 
magnetization of the I spins without affecting the R and S responses. The quiet version 
of the soft 7r pulses was used (RE-BURP-Q) with a bandwidth of 40 Hz. The duration of 
the sequence (8r) was incremented to give phase steps of tt/2 radians over a range of 2tt. 

Fig. 4. (a) Experimental test of an 02 module which rotates the / magnetization through 
7r radians about an axis tilted through ir/4 radians in the xz plane, leaving the S and R 
spins largely unaffected. The RE-BURP-Q soft pulse had a bandwidth of 25 Hz and a 
duration of 220 ms. The soft pulse sequence of Fig. 1 was used, but with the assignment 
of soft pulses to S and R interchanged. In practice the tilt angle was made appreciably- 
larger than 7r/4 to compensate for spin-lattice relaxation effects, (b) Simulated spectrum 
allowing for strong coupling but neglecting relaxation. 

Fig. 5. Experimental test of the 03 module in which Jjs causes phase divergence of the I 
and S transverse magnetization spanning a range of ±2ir radians in steps of ±7r/2 radians. 
RE-BURP pulses with a 40 Hz bandwidth and a duration of 121 ms were used. In traces 
(b) and (d) the receiver phase was shifted by tt/2 for the / and S responses in order to 
reset them to pure absorption. 
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